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^ SUMMARY 

The  wind  tunnel  wall  interference  in  transonic  speed  is  formulated 
as  a perturbation  to  the  basic  flow  around  the  airfoil  in  free  air.  The  pertur- 
bation equation  is  derived  from  the  transonic  small  disturbance  equation 
and  is  linear  but  with  variable  coefficients  containing  the  non-linear  solution 
of  the  basic  flow.  The  equation  is  solved  numerically  by  a direct  matrix 
method  using  the  classical  boundary  condition  for  a porous  wall.  The 
solution  in  terms  of  lift  versus  angle  of  attack  agrees  well  with  that  calcu- 
lated directly  from  the  small  disturbance  equation. 

\ 


RESUME 

L ’interference  que  causent  les  parois  d’une  soufflerie  transsonique 
est  exprimee  comme  une  perturbation  agissant  sur  Pecoulement  autour  de 
la  surface  portante  dans  Pair  libre.  L ’equation  de  perturbation  est  etablie 
a partir  de  l’equation  des  petites  perturbations  aux  vitesses  transsoniques; 
elle  est  lineaire  mais  comprend  des  coefficients  variables  contenant  la  solu- 
tion non  lineaire  de  l’ecoulement  primaire.  L’equation  se  resout  numeri- 
quement  au  moyen  d’une  methode  matricielle  directe,  en  utilisant  la  con- 
dition limite  classique  pour  une  paroi  perforee.  La  solution  exprimee 
en  termes  de  la  variation  de  portance  en  fonction  de  Tangle  d’attaque  Con- 
corde bien  avec  celle  calculee  directement  a partir  de  1 ’equation  de  faible 
perturbation. 

Aceessi  . 

ntTT 

DD2  t 


CONTENTS 


Page 

SUMMARY  (hi) 

1.  INTRODUCTION  1 

2.  GOVERNING  EQUATIONS  AND  CONDITIONS  2 

3.  PERTURBATION  EQUATIONS  AND  CONDITIONS  3 

4.  INTERFERENCE  CORRECTIONS  4 

5.  METHOD  OF  SOLUTION  5 

6.  RESULTS  AND  DISCUSSION  8 

7.  CONCLUSION  9 

8.  REFERENCES  9 

LIST  OF  SYMBOLS  11 

APPENDIX  Coefficients  of  Finite  Difference  Equation,  Equation  (23)  21 

LIST  OF  FIGURES 

Figure  Page 

1 Problem  Layout  and  Co-ordinates  13 

2 Comparison  of  Different  Results  for  CL  vs  a for  NACA  0012 

Airfoil  in  Perforated  Wind  Tunnel  and  Free  Air  at  M = 0.80, 

H*/c  = 2 14 

3 Comparison  of  Pressure  Distributions  at  the  Wind  Tunnel  Walls 

for  NACA  0012  Airfoil  at  M = 0.8,  H*/c  = 2,  P*  = 0.47  and 
o = 1°  (Geometrical)  15 

4 Comparison  of  Free  Air  Pressure  Distributions  at  y = ±H  for 

NACA  0012  Airfoil  at  M = 0.8,  or  = 1° 16 

5 Corrections  for  Angle  of  Attack  and  Mach  Number  from  Present 

Method  and  the  Linear  Theory  17 

6 Pressure  Distributions  at  the  Wind  Tunnel  Walls  for  a Transonic  Airfoil 

at  M = 0.6,  a = 1.63°,  H*/c  = 3,  P*  = 1.50,  P*  = 0.50  18 

7 Pressure  Distribution  at  the  Wind  Tunnel  Walls  for  a Transonic  Airfoil 

at  M = 0.6,  a = 3.90°,  H*/c  = 3,  P*  = 1.50,  P*  = 0.50 

(iv) 


19 


.1 


A PERTURBATION  THEORY  OF  TWO-DIMENSIONAL  TRANSONIC 
WIND  TUNNEL  WALL  INTERFERENCE 


1.  INTRODUCTION 

Interference  of  the  wind  tunnel  wall  on  the  flow  around  a model  is  an  important  and 
practical  subject  in  aerodynamics  of  simulation.  Through  the  years  a large  amount  of  effort  has  been 
devoted  to  the  problem.  An  excellent  summary  of  the  results  of  these  efforts  can  be  found  in  the 
review  by  Garner  et.  al.  (Ref.  1).  The  classical  approach  to  the  problem  is  to  consider  the  wall  inter- 
ference as  a perturbation  to  the  basic  flow  past  the  model  in  free  air.  The  flow  induced  by  the  perturba- 
tion is  then  treated  as  corrections  to  the  measured  results  of  the  model.  This  method  works  well  for 
subsonic  flow  because  the  linear  nature  of  the  flow  regime  warrants  the  principle  of  superposition  of 
solutions. 


In  the  transonic  range,  the  governing  equation  for  the  flow  past  a model  is  non-linear.  The 
classical  method,  in  principle,  can  no  longer  be  applied  and  elaborate  calculations  of  the  flews  bounded 
by  tunnel  walls  because  necessary  (Refs.  2,  3,  4).  These  rigorous  methods  usually  involved  numerical 
solutions,  require  a large  amount  of  computation  and  thus  are  not  favourable  for  practical  applications. 
The  subsonic  theory  has  also  been  applied  in  the  transonic  range  (Ref.  5)  based  on  the  fact  that  the  far 
field  structures  for  the  subsonic  and  the  transonic  flows  are  highly  similar.  However,  whether  the 
subsonic  linear  solution  can  be  applied  all  the  way  to  the  airfoil  around  which  the  flow  is  definitely 
non-linear,  has  not  been  proven. 

The  simplicity  of  the  classical  method,  on  the  other  hand,  is  attractive  to  practical  engineering. 
It  would  be  desirable  if  a method  could  be  derived  with  the  preservation  of  the  simplicity  of  the  classical 
method  and  at  the  same  time  being  correctly  formulated  for  the  transonic  range.  The  present  report  is 
a first  attempt  in  this  direction. 

In  keeping  the  simplicity  of  the  classical  method,  the  present  approach  considers  that  the  wall 
interference  is  a small  perturbation  to  the  basic  flow  around  the  airfoil.  However,  the  basic  flow  is  now 
transonic  and  non-linear.  Within  the  framework  of  transonic  small  disturbance  theory,  the  perturbation 
equation  is  then  derived  from  the  transonic  small  disturbance  equation.  The  resulting  equation  is  now 
linear  because  of  the  small  perturbation  assumption  but  with  variable  coefficients  containing  the  non- 
linear solution  of  the  basic  flow.  This  equation  is  then  solved,  satisfying  the  imposed  boundary  condi- 
tions at  the  tunnel  wall  and  the  airfoil.  The  solution  thus  gives  the  first  order  correction  to  the  basic 
flow  due  to  the  presence  of  wind  tunnel  walls.  The  results  of  the  present  method  compare  favourably 
with  those  directly  calculated  from  the  small  disturbance  equation  (Refs.  2,  4).  This  comes  as  no 
surprise  since  the  success  of  applying  perturbation  techniques  to  the  transonic  non-linear  basic  flow  has 
been  demonstrated  by  Cheng  and  Hafez  (Ref.  6)  for  three-dimensional  flows  and  Nixon  (Ref.  7)  for 
two  dimensional  cases. 

The  present  method  is  simple  to  use  for  practical  wall  interference  calculations  as  the  linear 
equation  can  be  solved  by  a standard  routine  with  a short  computing  time.  This  is  particularly  important 
if  parameters  such  as  tunnel  wall  porosities  have  to  be  determined  by  matching  the  solution  with  the 
measured  wall  pressure  distributions  (Ref.  5).  For  this  case,  the  non-linear  equation  of  the  free  air  basic 
flow  will  be  solved  only  once,  and  the  searching  routine  involves  solving  the  linear  equation  only. 

The  present  calculation  also  shows  that  the  angle  of  attack  correction  obtained  by  the 
subsonic  linear  theory  applied  in  the  transonic  range,  is  surprisingly  close  to  that  of  the  present  results. 
The  Mach  number  correction  from  the  linear  theory  however,  does  not  compare  so  favourably.  This  is 
attributed  to  the  fact  that  in  the  linear  approach,  the  non-linear  contribution  to  the  doublet  strength, 
which  is  strongest  near  the  airfoil,  is  completely  neglected. 
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2.  GOVERNING  EQUATIONS  AND  CONDITIONS 

The  problem  considered  is  a transonic  inviscid  flow  past  an  airfoil  model  situated  in  a wind 
tunnel  (Fig.  1).  The  airfoil  has  a chord  of  unity  and  a thickness  5.  The  tunnel  walls  are  located  at  a 
distance  H above  and  below  the  model  (the  model  need  not  be  located  at  the  middle  of  the  tunnel; 
the  central  position  is  chosen  for  simplicity).  The  airfoil  is  assumed  to  be  thin  and  the  angle  of  attack 
small.  Within  the  framework  of  the  transonic  small  disturbance  theory  (Refs.  8,  2),  the  equation 
governing  the  flow  inside  the  tunnel  is  thus 

(K  - (7  + 1)  0X  ] 0XX  + 0yy  = 0 (1) 

K is  the  transonic  similarity  parameter  defined  as 


1 - M2 

M2m  52/3 


(2) 


where  6 is  the  thickness  of  the  airfoil  and  M the  free  stream  Mach  number.  The  velocity  potential  0 
and  the  co-ordinate  normal  to  the  flow  direction  y have  been  scaled  according  to  the  similarity  trans- 
formation. 


0 = 6’2/3  Mn  0* 
y = 6I/3  Mm  y* 

where  the  started  variables  denote  their  corresponding  quantities  in  the  physical  plane. 

The  boundary  condition  on  the  airfoil  is  that  the  flow  should  be  tangent  to  the  surface. 


(3) 


<0y)„  „ = M"'m 


u ,G 


<0F\  a 
Wu.e  " « 


at  y = ±0  on  airfoil 


(4) 


where  F is  the  profile  of  the  fairing  and  a the  angle  of  attack.  The  tunnel  wall  is  perforated  for  flow 
ventilation  as  required  for  the  transonic  tests.  The  boundary  condition  is  given  as  (Refs.  1,  2) 


0X  ± 0y  = 0 at  y = ±H 


(5) 


where  P is  the  porosity  of  the  tunnel  wall  and  H the  distance  between  the  airfoil  and  the  wall.  These 
parameters  are  also  seeded  as  follows: 


H = 61/3  Mm  H* 

P = 61/3  M'm  P* 

The  exponents  m and  n in  the  scaling  are  set  to  be  (Refs.  2,  9) 

m * 1/2,  n = 3/4 

If  there  is  shock  wave  occurring  in  the  flow  field,  the  weak  solution  should  satisfy  the  discontinuity 
condition  across  the  shock  (Refs.  8, 10) 
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<K-(7+1)0x>I<M2  + IIM2  = 0 (7) 

II 0 D = o 

where  the  shock  shape  is  expressed  as 

x = xD(y)  (8) 

The  symbol  |[  [)  denotes  the  difference  and  < > the  arithmetical  mean  across  the  shock. 

3.  PERTURBATION  EQUATIONS  AND  CONDITIONS 

In  the  present  formulation,  the  basic  flow  field  is  that  of  an  airfoil  in  the  free  air.  The 
restriction  of  the  tunnel  walls  is  considered  as  small  perturbation  for  the  basic  flow.  Thus  the  velocity 
potential  can  be  written  in  two  parts,  the  free  air  potential  and  the  interference  potential 

0 = 0 o + 01  <9) 


where  <pl  is  one  order  of  magnitude  higher  than  <pa . From  the  problem  defined  in  the  last  section,  the 
equations  and  the  conditions  for  the  free  air  and  the  interference  flows  can  be  derived. 

Free  air  flow 


[K-  (7+l)0ox] 

^oxx  +0 oyy  " 0 

(10a) 

(0oy)ue  = Mn'm 

(10b) 

0ox>0Oy  “*■  0 

as  X2  + y2  -*  oo 

(10c) 

and  <K  + (7+1)0ox 

> If  0OX  l2  + II  </>oy  D2  = 0 

(10d) 

o 

II 

o 

■o- 

(10e) 

Interference  flow 

[K-(7+l)0Ox 

l^lxx  - (7+  l)0oxx  0IX  +<t>lyy  = 0 

(11a) 

O 

II 

at  y = ±0  on  airfoil 

(lib) 

i 

0lx  * p 01y  “ 

-(«ox  * ^ </>oy)  at  y = ±H 

(11c) 

The  shock  condition  required  further  consideration  because  the  shock  location  may  be  shifted, 
shock  iocation  can  be  written  as 

The 

xD  = x0D(y)  + X!D(y) 

and  the  shock  conditions  should  be  satisfied  along  xD.  However,  calculation  of  the  first  order  equation, 
Equation  (11)  depends  on  the  zero  order  solution  and  its  shock  location.  Thus  the  shock  condition  for 
the  perturbation  equation  has  to  be  transferred  to  the  unperturbed  boundary  by  analytic  continuations 
(Ref.  10).  The  shock  conditions  for  the  interference  flow  is  therefore 
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<(7  + 1)  l0ix  (x0.y) + ^oxx  (xo.y)  xiD(y)l>  = 2 x^y  x?y 
ll^i  (xo.y)  + 0ox  (xo.y)  xiD(y)ll  = 0 


(12b) 


Far  Field  Boundary  Conditions 


In  computing  the  free  air  solution,  the  far  field  boundary  condition  has  to  be  specified.  If 
the  computed  flow  field  does  not  extend  to  infinity,  the  asymptotic  solution  of  the  transonic  small 
disturbance  equation,  Equation  (10)  can  be  used  (Refs.  2,  9) 


<MX- y) 


r K1/2y  D 

— tan"1  + 


2n 


2ttK1/2  x2  + y2 


(13) 


where  the  doublet  strength  D is  given  as 


1 ,0.5  . . 7 + 1 , ,+»  j 2 


D ■ 5/o:5,p"'f')<ix+  4 


/ / 0X  dxdy 


(14) 


and  f is  the  circulation  scaled  as 


r = 5 2/3  Mn  r* 


By  satisfying  the  Kutta  condition,  the  circulation  P is  equal  to  the  potential  difference  across  the  airfoil 
plane  at  the  trailing  edge.  These  conditions  should  be  incorporated  with  the  zero  order  equation  and 
the  other  conditions  given  in  Equation  (10)  to  complete  the  formulation  for  the  free  air  case. 


For  the  interference  potential,  the  conditions  at  far  upstream  and  downstream  have  to  be 
specified.  The  asymptotic  solutions  for  this  case  are  given  by  Murman  (Ref.  2)  and  Catherall  (Ref.  4). 

If  the  flow  field  extended  to  infinity  at  both  ends,  then  the  condition  can  be  written  up  to  the  first 
order  as 

0,  = -0O  at  x -*•  ±°°  (15) 

where  0O  is  given  by  Equation  (13)  at  the  limit  of  x approaching  infinity.  The  potential  jump  at  the  trailing 
edge  will  give  the  change  of  circulation  due  to  the  wall  interference 

AT  = 0i+  - 0]"  at  the  trailing  edge  (16) 

The  formulation  for  the  interference  potential  is  thus  complete. 


4.  INTERFERENCE  CORRECTIONS 


In  the  classical  method  of  the  tunnel  wall  interference  (Ref.  1)  the  interference  potential  is 
calculated  without  enforcing  the  boundary  condition  on  the  airfoil  surface  by  considering  that  the 
airfoil  shrinks  to  a point  singularity.  The  correction  is  obtained  by  evaluating  the  velocity  component 
at  the  location  of  the  singularity.  The  stream  wise  component  is  considered  as  the  change  of  the  free 
stream  velocity  (or  Mach  number)  and  the  normal  component  the  change  of  angle  of  attack.  These  can 
be  written  as  follows  (Ref.  1) 


Act  = b Mm'n  0ly  (18) 

This  approach  is  justified  if  the  tunnel  height  2H  is  much  larger  than  the  chord  of  the  model.  Thus  the 
perturbation  parameter  e is  of  the  form  1/H  and  that  e « 1.  Within  the  limit  of  the  first  order  of 
perturbation  in  e the  classical  approach  can  be  adopted  for  the  present  formulation  as  well.  For  the 
changes  of  the  lift  and  surface  pressure  distribution,  however,  the  details  of  the  flow  field  around  the 
airfoil  are  required.  These  can  be  calculated  by  satisfying  the  boundary  condition  on  the  airfoil  surface. 
The  change  of  the  surface  pressure  can  then  be  evaluated  as 

ACp(x)  = S2/3  Mn  (-20lx ) at  y = ±0  (19) 

The  correction  of  lift  is  thus  from  Equation  (16) 

ACl  = 52/3  Mn  (2AO  (20) 

The  change  of  the  pitching  moment  can  be  calculated  from  the  surface  pressure  distribution,  (Eq.  (19)) 


ACm  = / ' AC  (x)  (x  + 0.25)  dx 


referring  to  the  quarter  chord. 


5.  METHOD  OF  SOLUTION 

The  zero  order  equation  and  the  conditions  governing  the  free  air  flow  have  been  extensively 
studied  and  computer  codes  have  been  developed  for  its  calculation  (Refs.  2,  9, 11).  For  the  cases 
calculated  in  this  report,  the  zero  order  solution  is  generated  by  a computer  code  available  in  the  Lab- 
oratory*. The  first  order  equation  can  be  solved  in  a similar  manner  as  the  zero  order  equation  by  a 
line  relaxation  method.  However,  since  the  first  order  equation  is  linear,  the  difference  equation  can  be 
solved  directly  by  an  existing  method  for  solution  of  a system  of  linear  algebraic  equations.  Some 
details  of  the  numerical  procedure  for  the  solution  of  the  first  order  equation  are  given  as  follows: 

Co-ordinate  Transformation 

The  x-co-ordinate  is  transformed  into  £ as 


such  that 


£ = 0 -►  ±1 


x = 0 -*■  ±°° 


where  s is  a parameter  for  scaling  x.  With  the  transformation  the  infinite  extent  of  the  x axis  is  brought 
within  a finite  domain  bounded  by  £ = ±1  and  y = ±H.  For  finite  difference  formulation,  the 
£-axis  is  divided  into  NX  points  and  y-axis  into  NY  points.  A finite  difference  equation  can  then  be 
written  for  each  grid  point  inside  the  domain. 


The  author  wishes  to  thank  D.J.  Jones  for  providing  the  zero  order  solutions  for  cases  calculated  in  this  report  and  for  stimulating 
discussion  on  the  course  of  the  program. 
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Finite  Difference  Equation 

In  deriving  the  finite  difference  equation  for  the  first  order  problem,  the  flow  field  is  divided 
into  the  subsonic,  the  supersonic,  the  parabolic  and  the  shock  regions  according  to  the  local  property 
of  the  zero  order  solution.  An  appropriate  fin:  difference  formula  is  assigned  to  each  region  accord- 
ingly (Ref.  11). 

The  types  of  flow  are  defined  as  follows: 


Subsonic 


Supersonic 


Parabolic 


Shock 


jO  — 

<t>XC  ~ 


0O  A.  ° 

i+l,j  ‘ #MJ  d£  K 


— < 


2A£  dx  7 + 1 


r°  .O 

,o  _ 0iJ  " *i-2J  d£  ^ K 
<Axb Z~T7 > 


2A£  dx  7 + 1 


0XC  > 


K 


7+1’ 


*xh  < 


K 


7+1 


Kc  < 


K 


7 + 1 ’ 


<b  > 


K 


7 + 1 


(22a) 


(22b) 


(22c) 


(22d) 


where  the  superscript  o denotes  the  zero  order  solution.  The  finite  difference  equation  can  now  be 
written  for  each  grid  point.  In  the  subsonic  region,  a central  difference  formula  is  used  because  of  the 
elliptical  nature  of  the  second  order  partial  differentia]  equation.  In  the  supersonic  region,  a backward 
difference  formula  is  employed  for  its  hyperbolic  nature.  At  the  parabolic  points,  the  operation  is 
achieved  by  setting 

K-  (7+  1)0OX  = 0 

and  at  the  shock  points,  the  operator  is  equivalent  to  the  sum  of  the  operators  of  the  subsonic  and  the 
supersonic  types.  These  operators  are  of  second  order  accuracy. 


With  these  operators,  the  interference  potential  equation,  Equation  (11a)  can  be  written  in 
a general  form  for  a grid  point  (i,  j)  as 


G 01-2  ,j  + A0Mj  + + c<t>  ij  + i,j+i  + = F (23) 

The  details  of  the  coefficients  for  each  type  of  the  flow  and  the  boundary  conditions  are  given  in  the 
Appendix.  For  all  grid  points,  Equation  (23)  forms  a system  of  linear  algebraic  equations  which  can 
be  written  in  a metric  form  as 


MX  = Y 


(24) 


where  M is  an  N X N matrix  of  block  tri-diagonal  form 
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M = 


Qn-iPn 

and  the  vector  X and  Y can  be  written  in  partitioned  form  as 


p,R, 


Q?  P2  P'2 


R 


N -1 


(25) 


i 


(26) 


These  matrices  can  be  further  identified  as 


C2,l  D2,l 

A2,l 

E2 ,1 

®2,2^2,2  • 

Ql  = 

• 

R1  = 

• 

• D2,NY-1 

E2  ,NY  ^2  ,N  Y 

A2,NY 

E2,NY 

and  so  on.  The  vector  X is  the  unknown  <p  and  the  vector  Y the  right-hand  sides  of  the  equation,  F. 
Numerical  Solution 

The  flow  field  bounded  by  the  wind  tunnel  walls,  the  upstream  and  the  downstream  infinities 
is  divided  into  a rectangular  grid,  with  the  following  specification. 

1,  . . . , 41;  A£  = 0.05 
Ay 

1,  ....  22;  — = 0.04S78 

H 

The  airfoil  is  located  at  the  middle  of  horizontal  lines  NY  = 11  and  12.  Afte:  the  solution  is  obtained, 
the  solution  at  the  airfoil  plane  is  calculated  by  extrapolation. 

The  values  of  the  first  order  solution  at  each  grid  point  are  obtained  by  interpolating  the 
known  solution  calculated  by  the  existing  code.  The  coefficients  for  the  difference  equation,  Equa- 
tion (23)  are  then  evaluated.  The  matrix  equation,  Equation  (24)  is  solved  by  a subroutine  LEQTIB 
in  the  International  Mathematical  and  Statistical  Libraries  (IMSL)  (Ref.  13). 


NX 
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The  calculated  examples  reported  below  were  performed  on  an  IBM  3032  computer.  The 
computing  time  is  about  5CPU  seconds  for  each  run. 


6.  RESULTS  AND  DISCUSSION 

The  present  formulation  of  the  interference  problem  bases  explicitly  on  the  assumption  of 
small  perturbation.  As  discussed  in  Section  4,  it  is  required  that  the  perturbation  parameter  e = 1/H 
be  much  smaller  than  unity.  However,  for  most  model-tests  in  a wind  tunnel,  the  value  of  the  parameter 
ranges  from  0.30  to  0.50  which  is  relatively  large  for  a first  order  theory.  In  principle,  higher  order 
terms  in  the  perturbation  series  should  be  included.  However,  by  comparison  with  the  results  of  the 
exact  small  disturbance  solutions  (with  tunnel  walls),  it  has  been  found  that  for  practical  purposes  the 
first  order  correction  is  adequate. 

The  test  case  chosen  for  the  present  method  has  the  following  conditions, 

Airfoil  - NACA  0012 

M = 0.8 

H*/c  = 2 

pu,L  = 0-47  (Ref.  14),  0.44  (Ref.  2) 

This  case  has  been  calculated  by  Murman  (Ref.  2)  and  Catherall  (Ref.  14)  using  the  small  disturbance 
theory.  These  results  on  the  lift  correction,  unfortunately  do  not  agree  with  each  other.  Details  of 
tunnel  wall  pressure  distribution  are  also  available  for  Catherall’s  solution  (Ref.  15). 

The  results  for  the  lift  corrections  are  shown  in  Figure  2.  The  free  air  solutions  for  the 
present  calculations  are  those  of  Jones  (Ref.  12)  which  are  in  good  agreement  with  Murman ’s  results. 

The  solution  for  the  wind  tunnel  case  are  obtained  by  adding  the  first  order  solutions  of  the  tunnel 
interference  to  the  free  air  solutions.  The  results  of  this  first  order  theory  agree  reasonably  well  with 
those  of  Murman ’s  direct  solutions,  though  the  perturbation  parameter  c/H  has  a value  of  0.5  which  is 
rather  large  for  the  first  order  theory.  Catherall’s  solutions  both  for  the  free  air  and  the  wind  tunnel 
are  higher  than  Murman ’s  and  the  present  results.  The  correction  of  Murman ’s  wind  tunnel  results  back 
to  the  free  air  using  the  subsonic  linear  theory  are  also  shown  in  the  Figure  for  comparison. 

The  pressure  distributions  on  the  tunnel  walls  are  shown  in  Figure  3 with  the  results  of 
Catherall  (Ref.  15)  and  Mokry  (Ref.  5)  for  comparison.  The  present  results  follow  that  of  Catherall’s 
closely.  The  linear  subsonic  theory,  however,  fails  to  predict  the  pressure  peaks  on  both  walls.  It 
might  be  argued  that  away  from  the  airfoil,  the  non-linear  effect  decays  rapidly  and  the  linear  theory 
should  represent  the  local  flow  reasonably  well.  This  indeed  is  true  if  the  tunnel  wall  is  located  far 
enough  away  from  the  model  so  that  the  asymptotic  solution  is  sufficiently  accurate  there  ( |y|  =s  6) 
(Ref.  9).  For  the  present  case,  the  location  of  the  wall  has  a value  of  0.88  in  y which  is  far  too  small 
for  the  asymptotic  solution  to  be  valid.  This  can  be  seen  by  comparing  the  pressure  fields  at  a distance 
H away  from  the  airfoil,  as  generated  by  the  transonic  small  disturbance  solution  and  the  asymptotic 
solution  of  the  linear  theory  for  the  airfoil  in  free  air  (Fig.  4). 


Since  the  computation  grid  is  chosen  for  the  calculation  of  the  interference  correction  as 
described  in  Section  4,  it  is  too  coarse  for  the  calculation  of  the  detailed  pressure  distributions  on  the 
airfoil  surface.  Further  refinement  of  the  numerical  calculation  is  needed  if  the  surface  condition  of 
the  airfoil  is  required. 


Two  more  cases  have  been  calculated  to  check  against  the  linear  theory.  In  one  case  the 
flow  is  completely  subsonic  and  the  other  case,  supercritical.  The  airfoil  has  a relatively  thick  transonic 
profile  and  the  conditions  for  these  two  cases  are  listed  as  follows: 
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16%  thick  transonic  profile 

0.599 

3 

1.50 

0.50 

1.63°  a = 3.90° 

0.282  CL  = 0.6 


The  porosity  factors  of  the  upper  and  lower  wall  are  chosen  to  give  the  best  fit  of  the  linear 
theory  to  the  experimental  pressure  distributions  at  the  walls  (Ref.  5).  The  calculated  results  are  shown 
in  Figure  5 with  those  of  the  linear  theory  for  comparison.  For  the  low  lift  case,  the  flow  field  is  com- 
pletely subsonic,  the  present  results  are  fairly  close  to  those  of  the  linear  theory  as  expected.  The  wall 
pressure  distributions  are  shown  in  Figure  6.  The  non-linear  compressibility  effect  can  be  observed  but 
is  not  too  pronounced.  At  higher  lift,  the  flow  over  the  upper  surface  is  supercritical  and  with  a shock 
terminating  the  supersonic  region.  The  present  results  give  a slightly  higher  (7% ) angle  of  attack  cor- 
rection and  much  larger  (25% ) correction  on  Mach  number.  The  wall  pressure  distributions  are  given 
in  Figure  7 in  which  the  effect  of  the  non-linear  compressibility  is  significant.  These  cases  again  show 
that  the  angle  of  attack  corrections  obtained  by  the  linear  theory  are  reasonably  close  to  the  non-linear 
results.  By  examining  the  far  field  asymptotic  solution  of  the  transonic  equation,  Equation  (13),  one 
may  observe  that  the  non-linear  term  appears  as  an  additional  doublet.  Thus  the  circulatory  flow  is 
very  close  to  that  of  the  linear  theory.  The  displacement  flow,  however,  will  be  affected  by  the  non- 
linear compressibility.  This  may  explain  that  the  Mach  number  correction  predicted  by  the  linear 
theory  has  a much  larger  discrepancy  than  that  of  the  angle  of  attack. 


7.  CONCLUSION 

The  transonic  wall  interference  is  considered  as  a perturbation  to  the  basic  flow  around  the 
airfoil  model  in  free  air.  Based  on  the  transonic  small  disturbance  theory,  the  perturbation  equation  is 
derived  from  the  non-linear  transonic  equation  and  is  linear  but  with  variable  coefficients  containing 
the  non-linear  solution  of  the  basic  flow.  With  the  boundary  conditions  imposed  at  the  tunnel  wall, 
the  equation  is  solved  numerically  by  a direct  matrix  method.  The  solutions  for  lift  versus  angle  of 
attach  agree  reasonably  well  with  those  directly  calculated  from  the  small  disturbance  equation.  The 
calculated  pressure  distributions  at  the  tunnel  walls  are  also  in  good  agreement  with  the  exact  results. 
The  present  method  is  convenient  to  use  for  practical  wall  interference  calculations  if  the  matching  of 
the  solution  with  the  tunnel  wall  pressure  data  is  required.  In  that  case  the  searching  routine  will  be 
solving  a linear  equation  only  which  is  much  quicker  and  simpler  than  solving  the  complete  non-linear 
equation. 


In  comparison  with  the  linear  subsonic  theory  which  is  commonly  used  for  correcting  wind 
tunnel  date,  the  present  results  show  that  the  angle  of  attack  and  the  Mach  number  corrections  obtained 
by  the  linear  theory  work  fairly  well  even  in  the  transonic  range.  The  linear  theory  underestimates  the 
corrections  by  a small  amount,  for  the  cases  here  presented,  which  is  well  within  normal  experimental 
accuracy. 
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LIST  OF  SYMBOLS 

CL  lift  coefficient 

Cp  pressure  coefficient 

D doublet  strength 

H tunnel  half  height 

K transonic  similarity  parameter 

M free  stream  Mach  number 

P porosity  factor 

s scale  for  x-axis 

x,y  Cartesian  co-ordinates 

a angle  of  attack 

7 ratio  of  specific  heats 

T circulation  around  airfoil 

5 maximum  thickness  of  airfoil 

£ transformed  co-ordinate 

0 velocity  potential 

I 


TUNNEL  WALL 


flow  direction 


TUNNEL  WALL 


FIG.  1:  PROBLEM  LAYOUT  AND  CO-ORDINATES 


lUM 


COMPARISON  OF  PRESSURE  DISTRIBUTIONS  AT  THE  WIND  TUNNEL  WALLS  FOR  NACA  0012  AIRFOIL 
AT  M = 0.8,  H*/c  = 2,  P*  = 0.47  and  a = 1°  (GEOMETRICAL} 


FIG.  5:  CORRECTIONS  FOR  ANGLE  OF  ATTACK  AND  MACH  NUMBER  FROM 
PRESENT  METHOD  AND  THE  LINEAR  THEORY 
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APPENDIX 

COEFFICIENTS  OF  THE  FINITE  DIFFERENCE  EQUATION,  EQUATION  (23) 


The  coefficients  of  the  finite  difference  equation,  Equation  (23)  are  derived  as  follows  for 
each  type  of  flow.  The  subscript  1 denoting  the  first  order  potential  has  been  dropped. 


Subsonic 


where 


with 


(Al) 


The  coefficients  A,  B,  C,  D and  E are  evaluated  at  the  grid  point  (i,  j)  and  the  functions  Z, , Z2  and  Z3 
are  introduced  because  of  the  co-ordinate  transform  Equation  (21).  A£  and  Ay  are  the  grid  sizes  along 
the  £ and  y direction  respectively. 
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At  the  wind  tunnel  walls,  the  boundary  conditions  Equation  (11c)  have  to  be  satisfied.  Along  the  lower 
wall,  at  grid  points  (i  j=l),  the  value  of  0;  j_}  is  replaced  by  the  boundary  condition.  Assuming  the  flow 
is  subsonic  there,  the  equation  becomes 


A0i.,,j  + C 0,j  + D0jj+,  + E0i+,  j = F at  y = -H,  j = 1 
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?L  is  the  porosity  factor  of  the  lower  wall. 


(A5) 


Along  the  upper  wall,  at  grid  points  (i  j-NY),  the  value  of  0;  j+  , is  eliminated  by  the  boundary 
condition.  The  equation  becomes 


where 


A0i.,  ,j  + B 0; j_i  + C0iJ  + E0i+l  j = F at  y = +H,  j = NY 


(A6) 


Pu  is  the  porosity  factor  of  the  upper  wall. 

At  upstream  and  downstream  infinity  of  the  flow,  the  boundary  condition  is  given  by  Equa- 
tion (15).  Thus  we  set 


0,  j = -0O  at  x - ±°°  or  i = 1,NX  (A7) 

j = 1,  2 NY 

In  the  classical  approach  of  tunnel  wall  interference,  the  interference  potential  i»  alculated  with  the 
airfoil  shrunk  to  a singularity  point.  The  boundary  condition  in  the  airfoil,  Equation  (1  lb)  need  not 
be  applied.  However,  if  the  details  of  the  flow  field  around  the  airfoil  are  required,  they  can  be  cal- 
culated with  the  boundary  condition  on  the  airfoil  enforced.  The  finite  difference  formula  of  the  second 
derivative  in  the  y-direction  applying  on  the  boundary  condition  can  be  written  as  (lief.  H) 


^yy)i,j  Ay2  [0i,j+l  0|,)] 


(0y ) 

’ y = » 


For  the  present  condition  in  Equation  (lib)  the  last  term  is  zero.  This  condition  should  be  applied  at 
both  the  supper  and  the  lower  surface  of  the  airfoil. 


